Photonic quantum memory in two-level ensembles based on modulating the refractive 
index in time: equivalence to gradient echo memory 
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We present a quantum memory protocol that allows to store light in ensembles of two-level atoms, 
e.g. rare-earth ions doped into a crystal, by modulating the refractive index of the host medium 
of the atoms linearly in time. We show that under certain conditions the resulting dynamics is 
equivalent to that underlying the gradient echo memory protocol, which relies on a spatial gradient 
of the atomic resonance frequencies. We discuss the prospects for an experimental implementation. 

PACS numbers: 



I. INTRODUCTION 

Quantum memory for light is an essential element for 
the photonic implementation of quantum communication 
and information processing [H, 0] • In recent years there 
has been a lot of work both on theoretical proposals and 
on experimental implementations [H, H| . Thus far opti- 
cal control, using relatively strong laser pulses, has been 
exploited for electromagnetically induced transparency 
and off-resonant Raman-type storage 0, Q in ensembles 
of three-level systems. More recently, a direct control 
of the transition dipole-moment has been proposed that 
emulates Raman-type quantum memories in a two-level 
atomic configuration 

In photon-echo based memories Q, the light- matter 
coupling is controlled in a more indirect way by exploiting 
the dephasing and rephasing of inhomogeneously broad- 
ened atomic ensembles. This includes the controlled 
reversible inhomogeneous broadening protocol the 
atomic frequency comb protocol 10j, and the gradient 
echo memory (GEM) protocol [11 1. The GEM protocol 
has allowed the demonstration of the highest memory 
efficiency (in the quantum regime) so far 12j. 

Recently Ref. [l3| proposed a quantum memory proto- 
col based on controlling the refractive index. Considering 
an ensemble of three-level atoms inside a host medium 
(e.g. rare-earth ions doped into a crystal), which is lo- 
cated in a circular optical cavity, they showed that a con- 
tinuous change of the refractive index of the host medium 
during an off-resonant Raman interaction between a sin- 
gle photon, a classical control pulse and the atomic en- 
semble allowed mapping the state of the single photon 
into a collective atomic excitation. 

Here we consider storing quantum states of light in an 
ensemble of two-level atoms in a host medium, where the 
refractive index of the medium can be modulated during 
the interaction of the light with the atoms. In contrast 
to Ref. [l3| here there is no optical control pulse (which 
is related to the fact that we consider two-level instead 
of three-level atoms) and no cavity. Interestingly, we find 
that under certain conditions the considered system leads 
to dynamics that are equivalent to those of the GEM pro- 



tocol; controlling the refractive index of the host medium 
in time can mimic the effect of the spatial frequency gra- 
dient present in GEM. 

This paper is organized as follows. In Sec. |TT]we derive 
the dynamical equations for our system under certain 
conditions, which we discuss in detail. In Sec. IIIII we 
compare these results to the GEM protocol. In Sec. HVl 
we discuss a possible experimental implementation of our 
protocol. Sec. V contains our conclusions. 



II. MAXWELL-BLOCH EQUATIONS 

Here we study the propagation of the light and its 
interaction with two-level atoms inside a host medium 
whose refractive index varies in time. We show that in a 
certain parameter regime the time-dependent refractive 
index does not play a role in the propagation equation 
for the light. In contrast, it plays an essential role in 
the dynamics of the atomic polarization. For simplicity, 
we assume the field is propagating in a certain direction 
with a fixed linear polarization. (This is well justified for 
our choice of possible implementation in a waveguide, see 
below.) The wave equation for the electric field operator 
is analogous to the classical equation, namely 
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where E is the electric field, z is the direction of propa- 
gation, fj,Q is the vacuum permeability, D is the electric 
displacement field, e is the permittivity of the propaga- 
tion medium and P is the polarization of the embedded 
two-level atoms. There are thus two fundamentally dif- 
ferent contributions to D. The eE term is due to the 
permittivity of the host medium, whereas P describes 
the polarization of the two-level atoms that are the ac- 
tual memory system for the light. 

Consider the case where e is time-dependent. The per- 
mittivity of the medium is related to its refractive index 
as e(t) — n 2 (t)tQ. We consider a medium with a linearly 
changing refractive index, n{t) — rii + id. Based on this, 
only the first derivative of the refractive index remains 
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in the Eq. ([T]), giving 



field can be written as 
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)E=—(2n 2 E + 4n(t)nE)+fi P, (2) 



H = 



where c is the speed of light. We now introduce the 
slowly varying components of the signal field E and the 
atomic polarization P, E = £ e -i("'ot-ko(*)*) and P = 

p e -t(u: t-k (t)z) _ Here the wave vectQr fc Q ( t ) =ki + kt = 

(rii + ht)^ 3 - is a function of time and k — S^a, where ojq 
is the central frequency of the signal. 

The wave equation can be greatly simplified provided 
that a number of (realistic) conditions are fulfilled (see 
also the appendix). The second-order spatial deriva- 
tive for £ can be dropped provided that the field ampli- 
tude changes appreciably over the length of the medium 
(such that the derivative is comparable to £ / L) and 
ko{t) 3> 1/L. Similarly, the second order time-derivative 
can be ignored if ujq 3> 1/r, where r is the duration 
of the pulse. The same conditions also allow one to 
drop the first and second order derivatives of the slowly 
varying polarization operator. We are interested in the 
regime where the extent of the pulse in space (outside 
the medium) , C = cr, is much greater than the length of 
the medium, L (nj Ej] ■ This allows one to drop the first- 
order time derivative of £ compared to the first-order spa- 
tial derivative. Finally, by assuming An <C n.i, where An 
is the total change in the refractive index, and kL <C 
one obtains the simplified propagation equation 
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This shows that, under the above conditions, the propa- 
gation equation remains unchanged compared to that of 
systems with constant refractive index (k does not play 
an appreciable role in the propagation), see e.g. fiol Il4|. 
The derivation of Eq. (3) from Eq. (2) is discussed in 
detail in the appendix. 

We now derive the dynamics of the polarization of 
the dopant atoms. The polarization of the j-th atom is 
pj = d\e>)<j3 gel where a ge — c/- 7 } (e 3 1 and {g 3 \d\e J ) is 
the matrix element of the corresponding dipole moment 
component between the ground and excited states. The 
collective atomic polarization at a certain position z is 
the sum over the individual atoms in a slice of width Az. 
The slow component of this collective operator is given 
by 



1 



i(u t— fco(t)zj) 



3=1 



, N 

(g\d\e)yCTge, 



where we assume equivalent dipole moment for all of the 
atoms; A and V are the cross-section area and volume of 
the light-atom interface; N is the number of the dopant 



atoms and a ge = J2j=i ^ge^ i s the average 

atomic polarization at position z. The Hamiltonian of the 
ensemble of the dopant atoms interacting with the light 
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where we assume uniform excited state energy Ml for all 
of the atoms. We can now derive the dynamics of the 
slowly varying collective atomic polarization using 
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For the next step we assume that all of the atoms are 
initialized in the ground state and the number of atoms 
N 3> 1. For weak (quantum) signals one can then ignore 
the change in the excited state population. 

Using Eqs. (|3I4|5I6|) . the above definition of the slowly 
varying field, and including the atomic excited state 
linewidth 7, one finds the Maxwell-Bloch equations de- 
scribing the interaction of the light with the collective 
atomic polarization in a medium with linearly time- 
varying refractive index, 



da ge (z,t) 
dt 



= -(7 + i(A + kz))a ge (z, t) + ig£(z, t),(7) 
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where £ 



and A = ft — cjq is the detuning. 



(3) The coupling constant g — (e\d\g) 



where uiq is 



the central frequency of the pulse and ei = n 2 eo is the 
initial refractive index of the medium. Note that the 
time dependence of the permittivity can be ignored in 
the definitions of £ and g because we are interested in 
the regime where An n^. 

The above set of equations shows the role of the lin- 
early changing refractive index of the host medium in the 
regime that we have discussed. One sees that the linear 
change of the refractive index in time results in a space- 
dependent frequency shift given by the kz term in Eq. 
([7]), see also Fig. QJb). The above Maxwell-Bloch equa- 
tions are identical to those underlying the GEM quantum 
memory protocol In the next section we therefore 

discuss in detail how the present proposal compares to 
GEM. 



III. COMPARISON WITH GRADIENT ECHO 
MEMORY 

The dynamical equations (8) derived in the previous 
section (under a number of realistic conditions) are ex- 
actly equivalent to those underlying the GEM quantum 
memory protocol pjj . In the GEM protocol, an initially 
narrow atomic absorption line is broadened by applying 
an external (longitudinal) field gradient. This longitudi- 
nal broadening allows one to accommodate all frequency 
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FIG. 1: (a) In the GEM protocol, a longitudinal energy shift 
in the atoms (solid dots) allows one to cover all of the fre- 
quency components of the incoming light, (b) In the protocol 
proposed here, due to the linear change of the refractive index 
in time, the light experiences an effective position-dependent 
frequency shift kz. This allows different frequency compo- 
nents of the light to interact on resonance with a spectrally 
narrow line of atoms. 



components of the incoming pulse, see Fig. HJa). Once 
the pulse is absorbed, the produced collective atomic ex- 
citation starts to dephase due to the position-dependent 
detuning. By inverting the external field that is used 
to generate the gradient, one can rephase the collective 
atomic excitation, which leads to the re-emission of the 
light. 

At first sight it may seem surprising that a time vari- 
ation of the refractive index leads to a spatial gradi- 
ent in the detuning. This can be understood in the 
following way. In the definition for the slow-varying 
field £, E = £ e -i{"ot-ko(t» , the fast- varying phase 
—i(u>ot — ko(t)z) has a temporal and a spatial part. One 
can define an effective frequency for the light by taking 
the time derivative of the phase, w e // = — kz. One 
can see that this corresponds to a spatial gradient in the 
frequency of the light, leading to a spatial gradient in the 
light-atom detuning, see Fig. [Tfb) . A quantum memory 
can then be realized in analogy with GEM. The light is 
absorbed while the refractive index is changing linearly in 
time. Then the refractive index is kept constant for stor- 
age. The light can be retrieved by changing the refractive 
index linearly in time again, but with the opposite sign to 
before. As for other two-level memory protocols includ- 
ing GEM, the storage time can be increased if necessary 
by transferring the population (after the absorption has 
been completed) from the excited state to a long-lived 
third state. 

The total shift in the effective frequency over the length 
of the medium is kL. In order to accommodate all fre- 
quency components of the signal kL has to be larger than 
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FIG. 2: The efficiency of the proposed memory protocol 
based on refractive index modulation in terms of the ini- 
tial optical depth di n . The efficiency is given by e~ 2jT (l — 
exp (— din'y /kL)) 2 . The figure shows the efficiency for differ- 
ent pulse durations r, relative to the excited state line width, 
7. We assume kLr = 2. Depending on the available optical 
depth, one can optimize the achievable efficiency by choosing 
an appropriate pulse duration. 

the frequency bandwidth of the signal, Auj = 1/r. There- 
fore, kL can be understood as the memory bandwidth. 
This is equivalent to the role of r]L in GEM, see Fig. Q] 
and [3. 

The efficiency of the present quantum memory pro- 
posal can be found by analogy to GEM, see Eqs. (3, 4) 
in [14|. Converting the equations of motion to the fre- 
quency domain, one finds that the transmitted pulse is 
attenuated by a factor of exp(— /3ir), where (3 — n . 
This implies that the optical depth of the system is 
d = 2/3tt = = d in -?-. Here d is the optical 

ck kL 

depth that is associated with the effectively broadened 
line, with the initial optical depth d% n = 27r nNg L . The 
retrieval efficiency is then given by (1 — exp (— d)) 2 e~ 2 ~ <T , 
see Fig. [2j Here we have assumed that the decay of the 
excited state only has an effect during absorption and 
retrieval, but not during storage. As mentioned before, 
this can be achieved e.g. by transferring the excitation 
to a third, longer-lived state. Hyperfine ground states 
in rare-earth doped crystals can have coherence times of 
many seconds [15j . 



IV. POSSIBLE IMPLEMENTATION 

We now discuss a potential experimental implementa- 
tion of our proposed protocol. We propose to use thulium 
ions doped into a lithium niobate waveguide. This sys- 
tem was used in a recent implement at ion of an atomic 
frequency comb memory [16| . Lithium niobate is an 
attractive host for the present proposal because of its 
electro-optic properties, see below. The thulium ions in- 
teract with near-infrared light at a wavelength of 795 
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nm. The transition is naturally inhomogeneously broad- 
ened. We propose to prepare an initial atomic lincwidth 
of 7 = 10 MHz by optical pumping, which is very realis- 
tic. We consider the case where kLr = 2. This assures 
that the memory bandwidth is large enough that it can 
accommodate the incoming pulse. Assuming L = 1 cm 
this leads to the requirement An w 5 * 10~ 5 . We choose 
the pulse duration r = I/IO7. The above values assure 
that w > 1/t, L < C = ct, kL < 2%, An < n { and 
ko{t) 3> as required for the derivation in section II. 
For the given parameter values one needs di n = 50 to 
achieve about 70% efficiency, see also Fig. [2] 

We consider the refractive-index modulation of the or- 
dinary optical axis of lithium niobate by a fast varying 
electric field. We consider the case where the crystal 
is clamped (spatially confined) and the temperature is 
a few Kelvin. Under these conditions the refractive in- 
dex of the ordinary axis n s» 2.26 [17|. The change in 
the refractive index through the electro-optical effect is 
governed by A(4r)# = J2k r ^kE k , where i,j = 1 is as- 
sociated with the refractive index change of the ordinary 
axis. This means that a time-dependent external field 
in a certain direction, E^, can impose a time-dependent 
refractive index for the ordinary axis if there exists a 
non-zero linear electro-optical coefficient for that direc- 
tion, nifc. For lithium niobate r U 3 ps 10 * 10~ 12 m/V 
and ni 2 * -3 * 10- 12 m/V [HI, EH- This leads to 
An w 6*10~ 5 under 0.85*10 6 -2.5*10 6 V/m electric field, 
depending on the direction of the field in the 2-3 plane. 
This is equivalent to applying 8.5 — 25V to a system that 
has 10/im thickness, comparable to the waveguide used 
in Ref. [Tsj ] . The maximum change in the refractive in- 
dex in lithium niobate is expected to be 10 -3 , which is 
limited by the breakdown electric field fl9| . 

Applying the external electric field to change the re- 
fractive index is potentially accompanied by level shifts, 
due to the linear Stark shift, for the atomic ground and 
excited states. On the other hand, for a certain type of 
dopant, by having the external electric field orthogonal to 
the difference between the permanent electric dipole mo- 
ment of the ground and excited states, one can keep the 
resonant frequency between these states unchanged. In 
our proposed system the permanent dipoles are aligned 
with the 3-axis [20J, |2l| , therefore the electric field should 
be applied along the 2-axis in order to avoid level shifts. 



V. CONCLUSION 

We have proposed a memory protocol based on vary- 
ing the refractive index of the host medium in time, and 
shown that in a certain regime it is equivalent to the 
GEM protocol, even though the latter is based on a spa- 
tial frequency gradient. One may wonder why no similar 
spatial gradient was seen in the protocol of Ref. fl3j . 
which also considered a time-varying refractive index, 
but in a Raman-type system, in contrast to the two-level 
ensemble considered by us. This can be understood by 



noting that in the scheme of Ref. [l3[ the refractive in- 
dex modulation causes a spatial gradient in the frequen- 
cies for both the signal and the control fields, such that 
the two-photon transition frequency remains unchanged. 
This holds for co-propagating signal and control. It might 
be interesting to consider counter-propagating signal and 
control fields in this context, for which no such cancela- 
tion would occur. This might result in a protocol similar 
to Raman-GEM 

We found that a relatively small modulation of the re- 
fractive index should be able to provide sufficient memory 
bandwidth. We proposed a potential implementation in 
lithium niobate waveguides doped with rare-earth ions. 
However, other implementations could also be consid- 
ered. In general, it may be easier to control the behavior 
of the refractive index in time, compared to the spatial 
control of the atomic transition frequencies required in 
standard GEM experiments. 
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VI. APPENDIX: WAVE EQUATION 
SIMPLIFICATION 



Here we explain in more detail how to obtain Eq. (3) 
from Eq. (2). Substituting E = £ e -K"ot-k {t)z) and 
P = -p e -i(uot-ko(t)z) ^0 Eq. |J3J) and canceling the fast 
varying phase gives 

(d 2 + 2ik (t)d z -k 2 (t))£ (8) 
( 9 2 _ 2i(wo _ kz)dt _ (wq _ ^) 2 )£ 

c l 

= ^{2h 2 - 4in(i)n(w - kz) + An{t)hd t )£ 
+u (d 2 - 2i(u Q - kz)d t - (u) - kzf)V. 

In the left hand side of the equation, the term feg(t)£ 

cancels " ^p ui^E. The second order spatial derivative of 
the field can be dropped compared to the first order, as- 
suming that the field changes appreciably over the length 
of the medium, L, and fco(t) 3> Similarly, the second 
order time derivative of the field and the first and second 
order time derivative of the polarization operator can be 
ignored compared to the first order derivatives if ujq 3> - . 
This simplification is valid as long as oj ^> kL. These 
conditions lead to 

n 2 (t) 

(2ik (t)d x H y-(2iu d t -2oj kz))£ (9) 

c l 

= 4-( 2 " 2 - Un(t)huQ + 4n(t)hd t )£ - «oWnP. 
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One can rewrite this equation and use ko(t) = n(t)uio/c 
to simplify it to 



(d z + ^-(d t +ikz))£ 
c 

in 2 2n 2in(t)h *Mo w o 

= c 2 fc (£) c 2 k (t) *' + 2k {t) 



(10) 



V. 



Provided that An <C n(t) (and using again ujq S> 1/t), 
one has « |gf 9 t £ « The condition 

An <C n(<) also allows one to drop ^£ compared to 
^^-d t £. This simplifies the wave equation to 



(d z + ^(d t + ikz))£ 

c 2ko(t) 



V. 



(11) 



Provided that the extent of the pulse in space, C = ct, 
is much larger than the length of the medium, L, one 
can ignore ^^-dt£ in comparison with d z £ , which is the 
dominant term for the slowly varying component of the 
field. Finally, assuming that kL <C ^f; the term m W fcz g 



also can be dropped compared to the dominant term, 
leading to 



Q Q_ WOp 

The coefficient of can be rewritten as, 



(12) 



2fc,(l + M) " 2fc, ^ fc,-^' 



(13) 



Under the condition An rii this can be well approxi- 

2 

mated by ^ fc " , which leads to 



d z £ 



V. 



(14) 



which is the standard wave equation under similar (real- 
istic) conditions in the absence of a time variation of the 
refractive index. 
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